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Abstract. We investigate the phonon-induced decoherence and dissipation in a donor-based charge quantum bit re- 
alized by the orbital states of an electron shared by two dopant ions which are implanted in a silicon host crystal. The 
dopant ions are taken from the group-V elements Bi, As, P, Sb. The excess electron is coupled to deformation potential 
acoustic phonons which dominate in the Si host. The particular geometry tailors a non-monotonous frequency distri- 
bution of the phonon modes. We determine the exact qubit dynamics under the influence of the phonons by employing 
the numerically exact quasi-adiabatic propagator path integral scheme thereby taking into account all bath-induced cor- 
relations. In particular, we have improved the scheme by completely eliminating the Trotter discretization error by a 
Hirsch-Fye extrapolation. By comparing the exact results to those of a Born-Markov approximation we find that the 
latter yields appropriate estimates for the decoherence and relaxation rates. However, noticeable quantitative corrections 
due to non-Markovian contributions appear. 

PACS. 3.67.L.X, 63.20.Kr, 03.65.Yz 



1 Introduction 

During the last decade it turned out that solid-state based nano 
structures are promising candidates for the realization of quan- 
tum information processing devices [1 1. The building blocks 
are quantum mechanical two-state systems (qubits) and some 
of the proposed designs have been realized in groundbreak- 
ing experiments, see Ref. 1 1] for a recent review on this field. 
Thereby, various approaches have been undertaken, ranging 
from superconducting flux and charge qubit devices to devices 
using the spin or the charge degrees of freedom of individual 
electrons in confined geometries. Aiming at an extreme minia- 
turization of solid-state devices down to the nm-scale, it has 
been proposed to implant individual dopant atoms in a semi- 
conductor crystal and to use nuclear spin states of buried phos- 
phorus dopants to realize a spin-qubit (Kane's proposal [2|). 
Complementary to the Kane architecture, the charge degree of 
freedom of a single electron shared by two donor atoms in a 
host crystal can be used for the coding of the logical infor- 
mation, as proposed in Ref. [3 4|. Thereby, the logical states 
|0) and |1) are realized by the charge states of the double- 
donor-system with the excess electron either located on the 
left or on the right donor. The transition between these states 
occurs via tunneling of the electron between the two dopants. 
The charge qubits can in principle be controlled efficiently by 
external electric fields, e.g., by an applied gate voltage. This 
property renders the proposed architecture attractive for real- 
izing control schemes with available fabrication and read-out 
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technologies |4|. Experimental progress for this kind of ion- 
implanted Si:P nanostructures has been reported recently Q. 

On the other hand, solid-state qubits suffer from the large 
number of degrees of freedom due to their embedding in a com- 
plex many-particle environment. The environmental decoher- 
ence and dissipation lead to a deterioration of the performance 
of quantum logic operations and also strongly influence entan- 
glement between qubits |6| necessary for quantum gate oper- 
ations. Various sources of decoherence include nuclear spins, 
phonons, and electromagnetic fluctuations in the host crystal. 
To gain a detailed understanding of the various decoherence 
mechanisms, realistic model calculations have to be performed 
which then allow to sort out the different contributions. In this 
work we concentrate on the influence of a phonon bath on 
the shared electron. To be definite, we consider a charge qubit 
formed by two group-V donors as proposed in Ref. |7|. One 
donor is formed by a phosphorus atom while the second donor 
will be one of the class {Bi, As, P, Sb}. The donor pair is as- 
sumed to be implanted in a silicon crystal host and share a com- 
mon electron. We consider linear acoustic phonons coupled to 
the electron and determine the dynamics of the charge oscil- 
lations between the two donors. Due to the particular geome- 
try, a tailored phonon environment is formed for the electron 
which depends non-monotonously on the phonon frequencies. 
In order to provide accurate quantitative results on the decoher- 
ence and dissipation rates, we apply the numerically exact iter- 
ative quasi-adiabatic propagator path integral (QUAPI) scheme 
1 8 9 1 . In particular, we have improved the widely used method 
by providing a recipe to completely eliminate the Trotter dis- 
cretization error. This allows to obtain fully convergent exact 
results by extrapolation to a vanishing Trotter increment 1 10 1. 
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An appealing alternative to extensive numerical studies are ap- 
proximative calculations which, for instance, rely on the weak 
coupling between the qubit and the environment. The most fa- 
miliar Born-Markov or weak-coupling approximation (WCA) 
II II yields to simple closed expression for the decoherence 
and relaxation rates. However, they apply for typical situations 
when the bath has a smooth frequency distribution 11 II . In our 
case, the environment is particularly shaped by the geometry 
leading to a non-monotonous bath spectral density. Hence, it 
is not a priori clear whether the widely used WCA is appro- 
priate and a careful check is desirable. By comparing the exact 
numerical QUAPI results with the approximate WCA results 
below, we will show that for realistic parameters, the WCA typ- 
ically yields the correct order of magnitude for the decoherence 
and relaxation rates. However, differences are noticeable when 
a quantitative comparison is made. We furthermore note that 
the calculated phonon decoherence and relaxation rates com- 
prise a fundamental upper limit for the coherence properties of 
this architecture which can hardly be overcome. 

The presented set-up is related to a double-quantum dot 
charge qubit realized in a GaAs semiconductor 1 12 1, where the 
geometrical constraints induce charge qubit oscillations with 
noticeable non-Markovian corrections due to the particularly 
shaped phonon environment. While piezoelectric phonons dom- 
inate in GaAs, we have to consider here the dominating defor- 
mation potential electron-phonon coupling since the Si crystal 
displays inversion symmetry. 



2 The model 

To study the influence of the phonons on decoherence and dis- 
sipation, we assume that the charge qubit is isolated from any 
leads. It is formed by a pair of donor atoms embedded in a sili- 
con substrate, which share a single excess electron I3I4I . To be 
specific, we consider the situation of one donor being a phos- 
phorous atom while the second one is an individual donor atom 
X chosen from the group X G {Bi, As, P, Sb} 0. Then, the 
two logic states |0), |1) of the charge qubit are defined by the 
electron residing either at donor 1 or 2, respectively. The total 
Hamiltonian is given in terms of the standard spin-boson model 

tmrj 

H = H S + H B +H SB , (1) 

where Hs is the two-state Hamiltonian for the charge qubit, 
Hb models the phonon bath and H$ b includes the electron- 
phonon coupling. 



2.1 Model for the charge qubit 

We represent the Hamiltonian of the charge qubit in the ba- 
sis of the two localized electronic states denoted as \L) = |0) 
and \R) = |1), each being the Is orbital of the left/right donor 
atom, the latter being placed at the origin and at the position 
de y , see Fig. [0 for a sketch of the geometry. The localized or- 
bital belonging to the right (left) donor is fully described by the 
position vector of the electron, i.e., = r and — r + de y , 
respectively. In addition, we allow for an external constant en- 
ergy bias e which for instance could be due to a nearby capac- 
itive gate. In terms of the Pauli spin matrices er^, the two-state 



Hamiltonian then reads 



Hs — hAa x + hea z 



(2) 



The two lowest lying energy eigenstates |i? ± ) are given as an 
(anti-)symmetric superposition of the localized states \L) and 
\R) such that = (\L) ± \R))/V2 with energies E± = 

=FZ\/2. The tunneling amplitude then follows as A = E + — E- 
and is a function of the donor distance d. 

In order to determine the tunneling amplitude A, we have 
to calculate approximate eigenvalues of the lowest symmet- 
ric and antisymmetric energy-eigenstate. In principle, rather 
highly elaborated methods are available for their calculation, 
including the anisotropic conduction band dispersion of sili- 
con, the valley orbit interaction and valley interference effects 
171 141 . The latter leads to an oscillatory behavior of the tun- 
neling amplitude A for increasing the donor distance d. No- 
ticeably, the oscillations are weak if the two donors are placed 
in the [100]-plane of the Si host |7|. However, we aim for a 
detailed and quantitative understanding of the electron-phonon 
decoherence mechanism and thus resort to the simplest straight- 
forward procedure to determine the tunneling amplitude which 
is the well established linear combination of atomic orbitals 
(LCAO) 831151 . This tight-binding method is very successful 
for determining the molecular orbitals for the H% -molecule 
but can easily be generalized to our model by introducing an 
effective Bohr radius 1 16 1. When we neglect the conduction- 
band anisotropy, we can assume that the localized states |£) 
(£ = L, R) are represented by the Is orbital of each donor. 
They read 



10 = J^se- re/at 



(3) 



where is the effective Bohr-radius of the donor £ |16| and 
= |rj |. In the following, the left donor is assumed to be the 
phosphorous atom, whereas the right donor is taken from the 
group-V donors {Bi, As, P, Sb}. Hence, we introduce the ratio 
p such that cir —pa^- 

To calculate the energy levels an ansatz for the wave func- 
tion for the (anti-)symmetric part is made and the overlap 
between the two wave functions is calculated, yielding the en- 
ergies for the (anti-)symmetric state. If energies are scaled in 
atomic units, they read 1 151 



E± = E[ ±) (d)+Ei ± \d). 



(4) 



Here, E[ ± \d) is the kinetic energy and E% \d) is the poten- 
tial energy, both being functions of the (dimensionless) donor 
distance d. They read 



Ef\d) 
4 ±} (d) 



l±e~ d (l + d- d 2 /3) 
l±e~ d (l + d + d 2 /3) 

l±2e- d (l + d) + (l/d) - (l/d+l)e- 2d 



l±e~ d (l + d + d 2 /3) 



(5) 



Due to the fixed positions of the donors, there is no need to min- 
imize the energy with respect to the donor distance, in contrast 
to analogous calculations for the TJ^-molecule. According to 
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the LCAO calculations typical tunneling amplitudes for a dis- 
tance of d = 7.06 nm (which corresponds to a separation of the 
two dopants by n = 13 lattice sites) follow as A rj 16 meV. 
This is qualitatively consistent with the results obtained from a 
more refined approach taking into account interference effects 
in the Si band structure |j7). 



2.2 Coupling to linear acoustic phonons 

The phonon bath is due to the silicon host crystal and is mod- 
eled as usual in terms of the bosonic operators 6 q as 



(6) 



with the phonon dispersion relation w q . The electron-phonon 
interaction reads 117 181 



H SB = hY(a$N L + a?N R )(bl + 6_ q ) . 



q 



(7) 



Here, N% — 0, 1 is the number of the excess electrons on the 
donor £, respectively, and a q = A q e~ lq r «F^(q). The coupling 
constant A q depending on the wave vector q is specified below. 
Note that the phonons can propagate in all three dimensions, 
and the electron-phonon coupling is not isotropic in general 
1 19 1 . To take care of the charge distribution in each donor we 
define a form factor according to 



F 5 (q) = y> 



(8) 



where rt^(r) is the charge density of the donor £. The coupling 
Hamiltonian is rewritten in the form 1181 



(9) 



H S B = ~<J Z Y 5q0q + k-q) , 



with g q = [A q (Fi(q) — Fp(q))]. The charge density distri- 
bution then follows directly from Eq. as n^(r) = |(r|£)| 2 , 
which leads to the form factors i*i(q) = /z,(q) and F R (q) — 
f R {q)e-^ de y with / c (q) - 16/[4 + (q^) 2 ] 2 . 

In this work we focus on linear acoustic phonons with lin- 
ear dispersion relation w q = s]q|, s being the sound velocity 
for silicon (s w 9 x 10 3 m/s) |20|. Since the silicon crystal has 
an inversion center there is no piezoelectric coupling between 
electrons and phonons, wherefore the dominating coupling is 
due to the deformation potential. Thus, the coupling constant 
reads 

D I h 
-q\ L x t , (10) 



2p m Vu> q 



where D is the deformation constant for silicon (D 8.6 eV, 
see Ref. 1211 ). p m is the mass density of silicon (p m w 2.33 x 
10 3 kg m~ 3 , see Ref. 1 20 1) and V is the volume of the unit cell. 

All the properties of the phonon bath can be captured in the 
spectral density defined as 



GM = £ |<? q | 2 <5(c - Wq ) . 



(11) 



Using Eq. fllOi and the definition of the form factors and taking 
into account the geometry, the sum over q can be transformed 
into a continuous integral which can readily be carried out. One 
then obtains the spectral density 



G(w) 



MD 1 



ir 2 p m hs 5 



(12) 



x 

-2(4 



\~s aL 



—pa-L 
s 



where jo 1S the spherical Bessel function. The spectral density 
is sketched in the inset of Fig. [2] The low-frequency behav- 
ior is superohmic according to G(ui — > 0) oc ui 3 , while in the 
high-frequency limit, it decays algebraically as G(u — * oo) cx 
ui~ 5 . The crossover between these two limits occurs on a fre- 
quency scale uo c — s/ap = t^ 1 , where ap is the radius of 
the phosphorus donor (ap — 1.22 nm, see Ref. [16 1), yielding 
u c = 2.46 THz, which corresponds to an energy of 10.17 me V. 
As we will see below, typical tunneling amplitudes A are com- 
parable to this energy scale. Thus, the frequency distribution 
of the bath is no longer monotonous in the range of the rel- 
evant system energies. As common approximative analytical 
treatments 1311 1 1 of phonon-induced decoherence typically in- 
volve a smooth frequency distribution, it is not a priori clear 
whether their results are applicable to this situation. Moreover, 
the used Born-Markovian approximation which neglects bath- 
induced correlations might not describe properly the dynamics. 
This can be seen from the autocorrelation function 1 1 1 1 of the 
bath, i.e., 



L(t)=L R (t) + iL I (t) 



i J duG{uj) 



cotn — — cos Lot — % sin u>t 



,(13) 



which is shown in Fig.|2] The typical width of the correlation 
function is comparable to the time scale A^ 1 w lq~ x of the sys- 
tem dynamics. The Born-Markov approximation corresponds 
to replacing the strongly peaked real part L R (t) by a 5-function 
with the corresponding weight while the imaginary part Li(t) 
is often neglected. However, since the geometry tailors a spe- 
cific structured phonon environment for the charge qubit, it is 
not clear from the very beginning that the Markovian assump- 
tion is valid. It is the main purpose of this work to investigate 
this issue and compare exact real-time path integral simulations 
with approximate weak-coupling (Born-Markov) results. 



3 The improved QUAPI scheme 

The dynamics of the charge qubit is described in terms of the 
time evolution of the reduced density matrix pit) which is ob- 
tained after tracing out the bath degrees of freedom, hence 

p(t) =ti B K(t,0)W(0)K- 1 (t,0), 
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K(t,0) = Tcxp j - ^/ dt ' H 



K(t, 0) denotes the propagator of the full system plus bath, T 
is the time ordering operator and H is the Hamiltonian, Eq. Q. 
The full density operator W(Q) at initial time t — is as usual 
assumed factorizing according to W(0) cx p(0) exp[— Hb /(ks 
where ks is the Boltzmann constant and T is the bath tempera 
ture [ 1 1 1 . In this work the qubit dynamics always evolves from 
the initial state p(0) = \L) (L\. 

In order to investigate the dynamics of the system, we use 
the quasi-adiabatic propagator path integral (QUAPI) scheme 
1 8 1 being a numerically exact iteration scheme which has been 
successfully adopted to many problems of open quantum sys- 
tems 1 91121 . For details of the iterative scheme we refer to pre- 
vious works 1 8 9 1 and do not reiterate the scheme here again. 
However, we have improved the method at one important step 
and we will describe this in greater detail next. 

The algorithm is based on a symmetric Trotter splitting of 
the short-time propagator K(tk+i,tk) of the full system into 
a part depending on Hs and Hb + H$ b describing the time 
evolution on a time slice At. This is exact in the limit At — > 
but introduces a finite Trotter error to the propagation which 
has to be eliminated by choosing At small enough that conver- 
gence has been achieved. On the other side, the bath-induced 
correlations being non-local in time are included in the numer- 
ical scheme over a finite memory time r mem = KAt which 
roughly corresponds to the time range over which the bath au- 
tocorrelation function L(t) given in Eq. Jl 31 is significantly dif- 
ferent from zero. Note that for any finite temperature L(t) de- 
cays exponentially at long times [ 1 1 1 justifying this approach. 
To obtain convergence with respect to the memory time, K has 
to be increased until converged results have been found. How- 
ever, the numerical effort grows exponentially with the mem- 
ory length K and for the present two-level system, the memory 
length is restricted to typical values of K = 12 on a standard 
processor with 2 GB RAM for practical reasons. 

Thus, the two strategies to achieve convergence are coun- 
tercurrent. To solve this, the principal of least dependence has 
been invoked 1 9 1 to find an optimal time increment in between 
the two limits. However, here we show that the algorithm can 
be improved by applying a different strategy. 

We first choose some small enough time increment At. 
Then, one has to increase the memory time r mem by increas- 
ing K until convergence has been achieved. Typical results of 
this memory convergence check are shown in Fig. [3] Shown is 
the decoherence rate 7 for increasing memory time for differ- 
ent donor distances (p = 1) for the symmetric qubit e = 0. 
Note that the decay rate has been obtained by fitting the re- 
sults for the population difference P(t) — (<J z )t to an exponen- 
tially decaying cosine. The remaining error is the Trotter error. 
However, following Ref. 1101 . for any Hermitian observable, 
this symmetric Trotter error vanishes quadratically in the limit 
At — > 0. This opens the possibility to extrapolate the results 
to At — > 0, thereby completely eliminating the Trotter error. 
This is done by decreasing At from the initial value and then 
by finding the extrapolated exact result (of course, convergence 
has to be verified again for the smaller values of At). Typical 
results of this extrapolation procedure are shown in Fig.|4] in- 



dicating that the numerical values follow a line for decreasing 
(14) step sizes. Note that we consider P(tfi X ) at an arbitrary time 
tfix = 34.1a; c in this example. Indeed, we find the predicted 
behavior for the Trotter error to vanish and perform a linear re- 
gression to At — > 0, also shown in Fig.|4] The y-axis intersec- 
tion gives the numerical exact value for the observable of inter- 
est, in this case afflicted with a tiny error bar coming from the 
T)jjnear regression. In general, the convergence properties of an 
observable strongly depend on the involved parameters, similar 
to path-integral quantum Monte-Carlo simulations |22|. Differ- 
ent observables show different behaviors with decreasing Trot- 
ter step size At, as for instance the particle density in contrast 
to the energy of the system in Ref. 1221 . 



4 The dynamics of the charge qubit 

Equipped with the numerically exact improved QUAPI scheme, 
we can now study the dynamics of the charge qubit in detail. To 
extract the decoherence rate 7, the relaxation rate j r , the equi- 
librium population difference P^ and the oscillation frequency 
f2, we fit a combination of exponentially decaying cosine and 
sine functions 1 1 1 1 to the numerically exact data, from which 
the Trotter error has been eliminated. We can then investigate 
the dependence of the above quantities on the experimentally 
relevant parameters. We emphasize again that realistic assump- 
tions on the geometry of the system enter the spectral density 
Eq. (II 31 and thus allow to calculate quantitative realistic re- 
sults. 

One of the major goals of this work is to verify the Born- 
Markov (weak-coupling) approximation, since the later results 
in very simple and compact formulas for parameters govern- 
ing the dynamics. Hence, we compare the exact QUAPI results 
with results obtained within a WCA which are known as II II 



7 = 




(15) 


lr = 




(16) 


Q 2 = 


4A 2 eff [l - 2 Reu(2iA b )} + 4e 2 , 


(17) 


Poo = 


A b \ 2 J 


(18) 



The spectral function S(uj), related to the phonon spectral den- 
sity, Eq. (TO), via S(u) = G(ui) coth(hu/(2k B T)), repre- 
sents emission and absorption of a single phonon and A b = 
2 a/ A 2 s + e 2 is twice the effective qubit frequency. A c g is the 
effective tunnel matrix element at T = II II . which includes 
the renormalization by a Franck-Condon factor stemming from 
the high-frequency modes of the reservoir 1 11]. For our case, 
one easily finds that A c g « A with a deviation of less than 1%. 
The function u(z) is defined in terms of the frequency integral 



u{z) 



U! 2 + Z 2 



coth 



hie 

2k B T 



(19) 



J. Eckel, S. Weiss, and M. Thorwart: Phonon-induced decoherence and dissipation in donor-based charge qubits 



5 



4.1 Coherent charge oscillations for the symmetric 
qubit e = 

For the symmetric qubit with zero bias (i.e., only decoherence, 
no dissipation), we have calculated the time evolution of P(t) 
and have observed coherent charge oscillations. In order to quan- 
tify them, we define the quality factor Q = f2/(nj) where 
the frequency Q and the decoherence rate 7 have been been 
obtained from the fit as described above. We have performed 
extensive simulations for three different donor distances d for 
various combinations of donor atom species and show the re- 
sults as a function of the tunneling amplitude A in Figs.|5]|6| 
and0 each for a fixed donor distance d. A variation of A for 
a fixed donor distance can, for instance, be achieved by a small 
additional gate voltage which slightly distorts the Is orbitals 
leading to an increased overlap of the wave functions. 

For the smallest donor distance d = 4.34 nm, we observe in 
Fig.[5]that Q increases monotonously for increasing A. Thereby, 
the results for Q vary over two orders of magnitude for the dif- 
ferent donor species at large A. Moreover, the combination of 
two P donors or of one P and one (very similar) Sb donor dis- 
plays the best decoherence properties. The dashed lines in Fig. 
|5]display the results of the WCA given in Eqs. d!5l > and d!7l >. A 
reasonable agreement is found in this case. 

For intermediate donor distance d = 7.06 nm, see Fig.|6j 
Q first decreases but then increases again with increasing A. 
This can be understood by the fact that d determines the shape 
of the spectral density and, in particular, the location of the 
frequency of the cross-over, relative to the qubit frequency A. 
For the overall performance, the similar observation as for the 
smaller distance (see above) apply. Also in this case, the WCA 
seems to be appropriate although small deviations for all A can 
be observed which can be attributed to small non-Markovian 
corrections stemming from the specifically tailored phonon en- 
vironment. 

In the case of large donor distance d = 10.32 nm, see 
Fig. Q the differences between the various donor species al- 
most vanish and are only noticeable at small A. Also the WCA 
agrees well at large A and also yields the correct order of mag- 
nitude for small A although differences become noticeable in 
this regime. Note that in this case, Q decreases for increasing 
A, in contrast to the case of small and intermediate distances. 

Noticeably, we find that the Q-factor is independent of tem- 
perature for all relevant parameter combinations (not shown 
here). This is due to the fact that realistic temperatures cor- 
respond to frequencies of T = 6.5 x 10 9 Hz and hence all 
system frequencies are much larger. This behavior is in con- 
trast to what we have recently reported in GaAs DQD systems 
CD. 

Note that the oscillatory behavior of A for increasing d 1 7 1 
is not included in this simple LCAO approximation. However, 
when considering the Q~ factor in Figs.|51|6j and[7J the oscil- 
latory behavior of A for growing donor distances d does not 
affect Q substantially. This can be rationalized by considering 
the weak-coupling results Eq. ( 1151171 for e = 0. Then, it be- 
comes clear that the only part where A(d) appears is in the 
high-frequency part of G(uS) (assuming low temperature such 
the the coth approaches one and being interested in A » u> c ). 
The prefactors, which in principle contain A(d), drop out when 
the ratio is calculated. 



4.2 Dynamics of the biased charge qubit e^O 

When a finite bias e ^ is present, in addition to decoherence 
also relaxation occurs to a non-zero asymptotic value ^ 0. 
The corresponding decoherence and relaxation rates are also 
influenced by the presence of a bias in the sense that the effec- 
tive qubit frequency Aj, grows with increasing e. Then, the be- 
havior of the environmental frequency distribution is essential: 
if it grows with increasing frequency, decoherence and dissipa- 
tion will become more effective and if it decreases the environ- 
mental effects will diminish. This is what we observe from the 
results shown in Fig. [8] For comparison, we also show the cor- 
responding WCA results, which yield the qualitatively correct 
behavior while differences in the quantitative results occur. 



5 Conclusion 

To summarize, we have investigated the phonon-induced deco- 
herence and dissipation in donor-based charge qubits formed 
by a pair of donor atoms placed in a Si crystal host. The donor 
pair is formed by one P donor and one donor of the group 
Bi, As, P, Sb. We have employed the numerically exact quasi- 
adiabatic path-integral propagator in its iterative version. The 
major achievements of our work is twofold: (i) We have first 
improved the QUAPI scheme in the sense that the Trotter dis- 
cretization error can now be completely eliminated by extrapo- 
lating the results to vanishing Trotter increment, as it is known 
that the error vanishes quadratically. (ii) Beyond these method- 
ical aspects, we have obtained numerically exact results for 
the real-time dynamics of charge qubits under the influence of 
acoustic deformation potential phonons. Realistic assumptions 
on the tunneling amplitude enter via LCAO calculations of the 
wave functions and the qubit energies in our model. Moreover, 
we have included the particular phonon environment tailored 
by the particular geometry of the set-up via geometrical form 
factors and materials characteristics. No fitting parameters of 
any sort were utilized. 

In the absence of a static bias we have investigated the Q- 
factor of the charge oscillations as a function of the donor dis- 
tance and as well as a function of the tunneling amplitude. We 
have compared our results with those obtained from a WCA 
within an analytical approach in terms of real-time path-integrals 
and found that only small non-Markovian corrections appear. 
This can be attributed to the dominating super-Ohmic proper- 
ties of the phonon environment at small frequencies. Further- 
more we have investigated the dynamics in the case of a static 
bias and have found that the qualitative behavior of the deco- 
herence and damping rates follows the form of the environmen- 
tal frequency distribution. Non-Markovian corrections are also 
found in this case. 

At present, no experimental realizations of this setup is yet 
reported. Nevertheless, we emphasize that our results on the 
decoherence and dissipation induced by the electron-phonon 
coupling represent a fundamental upper limit to the coherence 
of such donor based charge qubits which can hardly be nego- 
tiated due to its intrinsic nature. This has to be seen in view of 
the DiVincenzo criteria 1 23 1 and also for the future realization 
of quantum information processes. However, the dominating 
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source of decoherence in this kind of qubit realization has to be 
investigated in realistic devices. 
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Fig. 1. Sketch of the geometry of a donor-based charge qubit formed 
two donor atoms at a distance d and the various angles of the phonon 
pro] 
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Fig. 2. The bath autocorrelation function (response function) L(t) = 
ifl(i) + iLi(t) for the spectral density G(u) (inset) of the phonon 
bath for the case of two P donors in a Si host (p = 1, deformation 
potential phonons), with s = 9 x 10 3 m/s, an effective Bohr radius of 
ap = 1.22 nm, and inter-donor distance d = 10.32 nm. The tempera- 
ture is T = 50 mK. 
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Fig. 3. (Color online) Check of convergence with respect to the mem- 
ory time Tmem ~ K At for the decoherence rate 7 (symmetric qubit 
e = 0) and for the donor distances d — 4.34 nm, d — 7.06 nm and 
d — 10.32 nm and the corresponding tunnel matrix elements obtained 
from the LCAO. The Trotter time increment is fixed to At = 0.55tJ c . 
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Fig. 5. (Color online) Quality factor as a function of the tunneling am- 
plitude A for different donor combinations and a small donor distance 
d — 4.34 nm. The symbols depict the exact QUAPI results while the 
dashed lines mark the results of the WCA. Temperature is fixed at 
T — 50 mK. 
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Fig. 4. (Color online) Example of the Trotter convergence for the the 
population difference of the qubit, P(tji x ), from which the quanti- 
ties of interest are extracted. In the lower sketch the tunnel amplitude 
was chosen as A/u) c = 2.24 and tf ix — 34.1o; c , and in the upper 
sketch A/ujc = 3.24 and tfi x = 18.2a;c. The memory-time is fixed 
to Tmem = 3.85/W and three values of K — 10, 11, 12 have been 
chosen. At T? nem /K 2 — > the value P(tfi x ) is shown as a result 
of the extrapolation At — > 0, with the error of the linear regression 
(horizontal bars). 
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Fig. 6. (Color online) Same as Fig. [3] but for an intermediate donor 
distance d = 7.06 nm. 
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Fig. 8. (Color online) Upper panel: Relaxation (7,-) and decoherence 
(7) rate for increasing bias e. Symbols are the exact QUAPI results 
while the dashed lines are the corresponding WCA results. Lower 
panel: Oscillation frequency O and asymptotic value P^. The re- 
maining parameters are d = 7.06 nm, A — 0.57oj c , and T = 50 
mK. 



